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Abstract. We report the first INTEGRAL timing and spectral analysis of the high mass X-ray binary source 
2S 0114+65 at high energies (5-100 keV). The pulse period was found at 2.668 hr with a high pulsed fraction, 
~80% in both the 20-40 keV and 40-80 keV energy bands. The spin-up trend over ~8 years was measured to 
be P = —8.9 x 10~ 7 s s _1 . The hard X-ray spectrum obtained with JEM-X/ISGRI is well described by a high 
energy exponential cut-off power law model where the estimated luminosity is 1.8 x 10 36 erg s _1 in the 5-100 keV 
energy band, for a source distance of 7.2 kpc. We tentatively identify a cyclotron resonance scattering feature at 
~22 keV with one harmonic, implying a magnetic field of ~ 2.5 x 10 12 G. 
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1. Introduction 

The X-ray source 2S 0114+65 is a high mass X-ray binary 
(HMXB) consisting of an accreting neutron star (NS) and 
a type Bla supergiant opti cal counterpart (L S 1+65010) 
at a distance of ^7.2 kpc l|Reig et al.lll996lh The pulse 
period was first measured at 2. 78 hr llFinlev et aDll992l) 
with a 11.59 day orbital period l)CramptoTi"et"lii"lll985jh A 
spin-up rate over ~11 yr was observed to be ~ 6.2 x 10~ 7 s 
s _1 (Hall et al. 2000). The long spin period makes this the 
slowest known X-ray pulsar. Different models have been 
suggested in the past to explain th e peculiar long period 



find evidence for a fluorescent iron line when the source is 
in a low emission state , consistent with measurements of 
lYamauchi et alJ ()l990|h 

In this paper, we report the results of high energy 
INTEGRAL observations of 2S 0114+65. The observa- 
tions are described in § [5] Imaging, timing, and spectral 
analyses are presented in § [3J A discussion of the results 
is given in § 0] 

2. Observations and Data Analysis 



suggested m tne past to explain trie peculiar long period 
pulsation seen in 2S 0114+65 (c.f. lHall et alJl200ot and 



Most recently, 



referen ces therein) 

l)l999(l suggest that such a slow period is possible if 2S 
0114+65 was born as a magnetar with an initial magnetic 
field strength of <T0 14 G, decaying to a current value of 
the order 10 12 G, allowing the NS to spin down to the 
measured spin period within the lifetime of its compan- 
ion. The magnetic field has not been measured up to now, 
as no cyclotron lines have been detected in the X-ra y spec- 
trum up to 20 keV. Recently, iFarrell et aD l)2005|) found 
evidence for a super-orbital period of 30.7 days, suggesting 
the presence of warped, precessing accretion disk. 

Observations by lHall et all l)200fj|) with RXTE show 
that the 3-20 keV spectrum is best fit by an absorbed 
power law with exponential high energy cut-off, with pho- 
ton index ~1.4, cut-off energy ~8 keV, and folding energy 
^20 keV, being typical parameters for a HMXB. They 



The present dataset was obtained duri ng the Target 
of Op portunity (ToO) A02 INTEGRAL l| Winkler et alJ 
Li k van den Heuvel| |2003^i observation of the Cas A region performed from 5-6 
and 7-9 December 2004 ( 53345.6 - 53346.8 and 53347.8 - 
53349.8 MJD), i.e. from part of INTEGRAL satellite rev- 
olutions 262 and 263. We use data from the coded mask 
imag er IBIS/ISGRI llUbertini et alJ l2003t iLebrun et all 
2003j) for a total exposure of 181.9 ks, and from the JEM- 



Send offprint requests to: E. W. Bonning 
* e-mail: erin.bonning@obspm.fr 



X monitor l|Lund et al.ll2003j) for a total exposure time of 
11 ks. For ISGRI, the data were extracted for all pointings 
with a source position offset < 12°, and for JEM-X with 
an offset < 3.5°. The spectrometer (SPI) was not used to 
extract the hard X-ray spectrum due to the lower sensi- 
tivity of this instrument with respect to IBIS /ISGRI for a 
weak source below 100 keV. Above 100 keV, 2S 0114+65 
was neither consistently detected in a single exposure nor 
in the total exposure time. Data reduction was performed 
usin g the standard Offline S cience Analysis (OSA) version 
4.2 l|Courvoisier et ai]l2003jh The background subtracted 
ISGRI light curve for the timing analysis was reduced us- 
ing dedicated software iiJight_extract (version to be in- 
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20-40 keV 



Fig. 1. The 20-40 keV IBIS/ISGRI mosaicked and de- 
convolved sky image of the ~ 182 ks observation. Image 
size is ~ 15° x 5°, centered at 2S 0114+65 position. The 
pixel size is 5'. 2S 0114+65 was detected at a significance 
of ~ 37 a. 

eluded in OSA 5.0). The ISGRI light curve per energy 
band was obtained for each 300 s long time bin. 



3. Results 

3.1. ISGRI Imaging 

Fig. ^ shows a significance map around the source 2S 
0114+65 in the 20-40 keV energy range. Single pointings 
were deconvolved and analyzed separately, and then com- 
bined in mosaic images. The source is clearly detected at a 
significance level of 36.6 a. In the energy range 40-80 keV, 
the significance level was 10.8 a. From the ISGRI data, 2S 
0114+65 was observed with the imaging procedure at the 
position aj2ooo = 01 h 18 m 02!55 and <5j 20 oo = 65°17'02'.'3. 
The source position offset with r espect to the catalog posi- 
tion l|Bradt fc McClintocklll983(l is 0.5'. This is within the 
90% confidenc e level assuming th e imaging source location 
error given bv lGros et al.l I 2003h . 



3.2. Timing Analysis 

The ISGRI 20-40 keV and 40-80 keV high energy light 
curves were extracted from the images using all avail- 
able pointings and are shown in Fig. [21 During the first 
^30 ks of the observation, the flux rises to a mean 
count rate in the 20-40 keV energy band of ^1.6 cts 
s -1 (^9.5x 10 _11 erg s _1 cm' 2 ). Similarly, but less signif- 
icantly, the flux in the 40-80 keV energy bands attains a 
count rate of ~0.5 cts s _1 (~2.6x 10 _11 erg s" 1 cm -2 ). The 
count rate is converted to flux assuming a power law with 
an exponential high-energy cut-off power-law model (see 
Table ^ Section I3.3|l . The hardness ratio of these energy 
bands as a function of time indicates that no significant 
spectral variability was detected. 

We searched for coherent pulsations of the source in 
the 20-40 keV energy band where we have the best statis- 
tics. Using the 300 s binned light curve of the source in 



40-80 keV 




'ttf 



Time since (53345 MJD) 

Fig. 2. INTEGRAL/lSGRl light curve of 2S 0114+65 
in the 20-40 and 40-80 keV energy bands. The ISGRI 
data are taken from the whole observation and have been 
converted to flux assuming a power law with an exponen- 
tial high-energy cut-off model (see Table Each point 
corresponds to a ~2.2 ks interval. 



the 20-40 keV energy band, we computed a Power Density 
Spectrum (PDS) in the frequency range between 6.5 x 10 -6 
and 1.6xl0 -3 Hz from Fast Fourier Transforms. In the re- 
sulting PDS, an evident signal is present at v — 1.04 x 10 -4 
Hz. The highest peak in the power spectrum corresponds 
to a nominal period of 2.67 hr. The accurate period was 
found using the phase-delay fitting method. We divided 
the data over 8 intervals of ~ 23 ks each and folded the 
light curve corresponding to each of the intervals using 
the nominal pulse period. The pulse p eriod was not sig- 
nificantly detected in the first interval. lHall et all ((2000) 
suggest that 2S 0114+65 is an eclipsing system, which 
may explain the very low source flux in this first inter- 
val. However, using their eph emeris, the orbital p eriod 
11.63 ± 0.007 days given bv ICorbet etaH l|l999li . and 
no orbital period evolution, we find the region of low 
flux in our observation to lie directly within the error 
of the expected location of the eclipse. For intervals 2- 
8, we found the pulse phase delays of each of the folded 
light curves with respect to the second interval. We fitted 
these phase delays with the first two terms of the equa- 
tion 56 = 



= 56 Q + Scj(t — t ) + icj(i — t ) 2 , where 56 is 
the measured pulse phase delay, t is the mid time of the 
observation, 56 is the phase delay at t , Slo is the de- 
viation from the nominal pulse frequency, and Co is the 
pulse frequency derivative. The best determination of the 
mean pulse period was found to be P = 2.668 ± 0.004 hr. 
Errors are at 1 a confidence level. The short observation 
time and low statistics did not allow us to determine or 
give an upper limit on the pulse derivative. 

The folded light curve (excluding the first interval 
where the pulse period was not detected above 3 a) is 
shown in Fig. [3J The pulse profiles show a single peak 
up to 80 keV simil ar to tha t obse rved previously at 
lower energies, c.f. lHall et all l)2000h and iGorbet et all 
jlQQflh . The pulsed fraction (PF), defined here as PF = 
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Table 1. Spectral best fit parameters. 



Pulse Phase 

Fig. 3. IBIS/ISGRI background subtracted light curve of 
2S 0114+65 folded at the best found NS spin period (2.668 
hr). The epoch To is arbitrary. From top to bottom the 
energy ranges are 20-40 keV and 40-80 keV, respectively. 
The pulse profile is repeated once for clarity. 
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Fig. 4. Spin history of 2S 0114+65. The f illed triangle 
is from this work. The o pen circle is f rom | Finlev et alJ 
<ll992h . fille d circles fromlCqrbet et alJ l|l999j) . and open 
square from lHall et"aD l)2000|k 



{Irnax ~ J rnin) / { T max + J rnin) , IS measured at 76 ±5% in 

the 20-40 keV energy range and 84 ±20% in the 40-80 
keV band. 



3.3. Spectral Analysis 

The spectral an alysis was done using XSPEC version 11.3 
l|Arnaudll996|) . combining the 20-100 keV ISGRI data 
with the simultaneous 5-20 keV JEM-X data. Due to the 
short exposure time of JEM-X, and therefore lower statis- 
tics, we rebinned this data in a 5 channel energy response 
matrix in order to constrain the ISGRI spectrum at higher 
energy. For ISGRI we used a 15 channel energy response 
matrix. A constant factor was included in the fit to take 
into account the uncertainty in the cross-calibration of 
the instruments. The factor was fixed at 1 for the ISGRI 
data. A systematic error of 2% was applied to the JEM- 



Dataset 


JEM-X/ISGRI 


Model 


PL X HlGHECUT 


r 




Scut (keV) 




Sfbid (keV) 


99 1+12.1 


X 2 /dof 


14.5/15 


Ls-ioo keV (erg s 1 ) 


1.8 xlO 36 



X/ISGRI spectra which corresponds to the current uncer- 
tainty in the response matrix. All spectral uncertainties in 
the results are given at a 90% confidence level for single 
parameters. 

The joint JEM-X/ISGRI (5-100 keV) spectrum was 
first fitted with a simple model consisting of a power law, 
which was found inadequate with % 2 /dof=49.5/17. In or- 
der to compare with previously reported measurements 
l)Hall et alJl200ff) . we fit a common model for HMXBs, 
namely a power law with exponential high energy cut- 
off. The photon flux density adopts the form f(E) = 
AE~ T 'e^ Ecut ~ ~ E " ' E! ° ld , where A is a normalization constant. 
This model gives a better fit with x 2 /dof=14.5/15, and the 
resulting best fit parameters T = 1.6, E cut = 9.0 keV and 
-Efold = 22.1 keV, are consistent within their error s with 
the previously published values of lHall et alll|2000|) . 

We note that there is a feature in the residuals to this 
fit (see Fig. [SJ around 40-50 keV. This is suggestive of a 
cyclotron resonance scattering feature (C RSF), which is 
obser ved in many accreting X-ray pulsars l|Coburn et alJ 
2002) at high energies. We attempted to test the hypoth- 
esis of a cyclotron feature in the spectrum by modifying 
the power-la w model with a two -harmonic cyclotron ab- 
sorption line l|Mihara et alJll99C|h This fit shows a CRSF 
at 22 keV with a first harmonic at 44 keV and line widths 
9.8 and 1.7 keV respectively. Although this model removes 
the dip feature in the residuals, the x 2 /dof of 5.4/10 does 
not improve the fit. 

Absorption by neutral hydrogen was not included in 
these models since previous reported val ues of the hydro- 



gen column density, - 3 x 10 22 cm" 2 (|Hall et alj feoOO) 



did not have effects above 5 keV. Similarly to the timing 
analysis, we excluded the first 14 pointings of the total 
ISGRI spectrum. The best fit parameters of the model 
together with their errors are reported in Table ^ The 
flux was converted to lu minosity assuming a distance of 
7.2 kpc ijReig et al.lll996f) . The fairly large errors reported 
for the best fit parameters E cut and the related Sfoid are 
substantially due to the low statistics of the JEM-X data. 
The unfolded spectrum and the residuals of the data are 
shown in Fig. 03 

4. Conclusions 

We have carried out a temporal and spectral analysis 
for the first time with INTEGRAL data of the HMXB 
2S 0114+65. We determined the spin period to be 2.668 
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Fig. 5. INTEGRAL/ JEM-X (5-20 keV) and ISGRI (20- 
100 keV) unfolded spectra of 2S 0114+65 with the best- 
fit power law and exponential high-energy cut-off model. 
Residuals between the data and model are shown in the 
bottom panel in units of sigma. 



hr, which indicates that the long-term spin-up of the NS 
is still continuing. In Fig. 0] we show the detailed pulse 
period history of the source from the first determina- 
tion in 1986 with EXOSAT. Since 1996, the spin period 
has decreased from ~2.73 hr to ^2.67 hr. The measured 
spin period derivative from MJD 50451 to MJD 53348 is 
P = —8.9 x 10 -7 s s " 1 which is ~l-4 times as large as 
the value reported by lHall et alJ (2000). The increase in 
P likely occurs from an increase of accreting matter, and 
shows that the spin-up is not a linear trend. The pulsed 
fraction is observed to be significant up to 80 keV, 76% in 
the range 20-40 keV and 84% in the range 40-80 keV. The 
detection of pulsed emission up to 80 keV and the mea- 
sured spin-up trend strongly support the conclusion that 
2S 0114+65 is a rotating neutron star where the pulsed 
photons arise from the polar cap on the NS surface. 

The energy spectrum of 2S 0114+65 from 5-100 keV 
is well fit by the standard model used for X-ray binary 
pulsars, i.e. a power law with exponential high energy cut- 
off. Best fit parameters give a photon index 1.6, cut-off 
energy 9.0 keV, and folding energy 22.1 keV. The X-ray 
luminosity calculated from the model flux is 1.8 x 10 36 erg 
s _1 from 5-100 keV, assuming the source distance to be 
7.2 kpc. 

If the presence of a cyclotron resonance scattering fea- 
ture with E cyc ~ 22 keV is confirmed, the magnetic field 
will be ~2.5xl0 12 G, assuming a 1.4 M NS with a ra- 
dius of 10 km and scattering in the accretion column re- 
gion above the polar cap on the NS surface. The magnetic 
field of the NS is given by E cyc = 11.6 B 12 (1 + z)- 1 
keV, where B\i is the magnetic field in units of 10 12 G, 
and (1 + z) is the gravitational redshift. This value is con- 
sistent for X-ray pulsars with cut-off energies around 10 
keV. Unfortunately, the statistics of the high energy data 
(with our exposure time) prevent us from confirming the 
existence of the CRSF. More extensive INTEGRAL ob- 



servations at high energy should allow this result to be 
tested with a longer exposure time. 
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